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Abstract 

A generalization of the Density Functional Theory is proposed. The theory developed leads 
to single-particle equations of motion with a quasi-local mean-field operator, which contains a 
quasi-particle position-dependent effective mass and a spin-orbit potential. The energy density 
functional is constructed using the Extended Thomas-Fermi approximation. Within the frame- 
work of this approach the ground-state properties of the doubly magic nuclei are considered. 
The calculations have been performed using the finite-range Gogny DIS force. The results are 
compared with the exact Hartree-Fock calculations. 



PACS numbers: 21.60Jz, 31.15Ew, 31.15Gy 



1 



1 INTRODUCTION 



The Hartree-Fock (HF) method is one of the most widely used approaches in nuclear physics. 
It is based on the concept of independent particle motion in the mean field produced by effec- 
tive nucleon-nucleon forces which are generally non-local and density dependent. The resulting 
equations of motion contain the non-local single-particle potential (SPP) which is determined self- 
consistently. The comprehensive study of nuclear ground state properties within the HF method 
has been performed with the zero-range Skyrme-like forces ||l], |2| . The status of this problem is not 
completely the same for finite range forces. The exact solution of the HF equations in this case is 
not an easy task, mainly due to the non-locality of the SPP. For example, the complete solution 
of the HF equations was carried out in Ref. for a finite range effective force with a Gaussian 
formfactor using a harmonic oscillator basis. On the other hand, the M3Y effective force with a 
Yukawa formfactor was employed in Ref. Q within the Campi and Bouyssy local approximation 
for the single-particle density matrix. 

Thus the problem of localization of the non-local SPP becomes actual. As it is known the 
non-local exchange Fock part of the SPP is determined using the non-local single-particle density 
matrix (DM in the following). If one approaches the DM in terms of only local quantities such 
as the particle density and kinetic-energy density, the corresponding HF exchange energy becomes 
a functional of these local quantities. The equations of motion obtained from the resulting local 
HF functional are second-order differential equations. It is important to note that they do not 
contain any integral operators which lead to difficulties in the general non-local case. For instance 
the Negele and Vautherin expansion of the DM |p and its modification by Campi and Bouyssy ||5| 
enables the HF energy to be expressed in the pointed functional form. Recently, another approach 
based on the Extended Thomas- Fermi method (ETF, see, for instance, Ref. has been proposed 
to this aim in Ref. 

An alternative, in a sense, approach to the mean field theory, which is widely used in 
applications to electron systems, is based on the Kohn-Sham (KS) Q method within the framework 
of the Density Functional Theory (DFT). The original version of this theory was developed in the 
pioneering paper of Hohenberg and Kohn (HK) where an energy functional depending only on 
the local particle density was considered (so we shall call it local DFT). Later on other versions of 
the DFT were proposed (see, for example, Refs. [pA| , p^). In particular, the non-local extension of 
this theory was discussed by Gilbert in Ref. where the functional dependence on the DM was 
included. 

The main merit of the KS scheme consists of the following: it gives a way to obtain 
single-particle equations of motion for the local DFT. These equations contain the local mean-field 
potential that has to be determined self-consistently. Notice that in contrast to the approximate HF 
method, the DFT yields, in principle, the exact ground-state energies (and the referred quantities) 
of the many-body system. Concerning the single-particle spectrum only the last occupied level 
has the exact physical meaning of the chemical potential in the DFT which is just the particle 
separation energy. 

There is one more important difference between the HF and KS methods that is revealed 
in the applications to nuclei. The radial-dependent effective mass and the spin-orbit potential are 
essential components of the HF approach in nuclear physics. These two quantities arise owing to 
the kinetic-energy density and spin density dependence of the HF energy functional. However, in 
the original KS method the effective mass is constant and equal to the physical mass and there 
is no spin-orbit potential because this method starts from a local energy density functional. At 
the same time at least the spin-orbit potential is necessary for the realistic description of nuclear 
properties. It is possible to introduce the kinetic-energy density and spin density dependence of 
the DFT energy functional formally. In this case, to derive the single-particle equations following 
the ideology of the KS method, one would assume that any kinetic-energy density and spin density 
can be produced by the many-particle wave function describing the non-interacting system in some 
external potential with a spin-orbit component. But, as opposed to the case of the local particle 
density (see Ref. 0]), this statement has not been proved. 

To include the radial-dependent effective mass and the spin-orbit potential in the consid- 
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eration in a rigorous way, one would have to use the non-local extension of the DFT and to derive 
single-particle equations of motion directly from the energy functional of this theory. However the 
straightforward application of the standard variational principle to the non-local energy functional 
leads to serious difficulties in view of the specific properties of the pseudo-Hamiltonian obtained 



(see Ref. [13| for details). 



The main goal of the present paper is to develop the modification of the non-local gener- 
alization of the DFT which would be free from the above-mentioned shortcomings of the non-local 
theory. To this aim we define an energy functional that depends on the DM produced by a deter- 
minant wave function (in what follows we call it by a Slater-determinant DM). Although this DM 
generally does not correspond to any interacting fermion system nevertheless, we will show that the 
minimum of this functional coincides with the exact ground-state energy of the interacting system 
under consideration. Applying the variational principle, we derive the single-particle equations of 
motion which, in contrast to the KS equations, contain a non-local SPP. This is described in the 
second part. 

In the third part the quasi-local reduction of the DFT is developed. Within the quasi- local 
DFT the energy functional depends on the local particle densities as well as on the uncorrelated 
kinetic-energy and spin densities. The single-particle equations, which are obtained by the mini- 
mization of this functional, contain the local SPP, the uncorrelated radial-dependent effective mass 
and the spin-orbit potential. 

In the fourth part we derive a semiclassical HF energy functional within the quasi-local 
scheme starting from the recently proposed expansion of the DM in the Extended Thomas Fermi 
method The explicit formulae for the energy functional, the SPP and the effective mass 
are obtained using the finite-range Gogny force |^, as described in the fifth part. The residual 
correlation term is taken phenomenologically. In the sixth section we apply our method to the 
description of the ground-state properties of some doubly-magic spherical nuclei. The main results 
are set out in the summary. In the Appendix 1 some auxiliary formulae for the SPP are given. 
In Appendix 2 we describe a simple method to take into account the two-body correction of the 
centre-of-mass motion to the binding energy. 



2 THE NON-LOCAL GENERALIZATION OF THE DFT 

Let us consider a system of interacting fermions. In the nuclear case we are interested in 
systems with two kinds of particles, namely neutrons and protons. Let H be the non-relativistic 
many-particle Hamiltonian. The explicit form of this operator is not important here. One can 
associate it with the usual formulae: 

where 

T = -y—Ai (2) 
^ 2m ^ ^ 



is the kinetic-energy operator, vf^^ is the bare nucleon-nucleon (NN) strong two-particle interaction, 
yCoui jg Coulomb force acting between protons and the dots note the many-particle interactions 
if needed. 

The HK energy functional [10|, which only depends on the local particle density n, can 



be defined within the framework of the constrained search method as follows (see, for example. 



13,11 



EHK[n] = inf < ^\H\^ > , (3) 

where > is an arbitrary normalized A^-particle state. The short notation ^ ^ n hereinafter 
means the many-to-one mapping of the wave function ^{xi, . . . ,xi\i) to the local density n(r) i.e. 
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it means that the following equalities are fulfilled 

n(r) = np{r) + n„(r) , (4) 

a 

p{x,x') = N '^{x,X2,...,XN)'^*{x',X2,...,XN)dX2---dXN , (6) 



where p{x,x') is the single-particle DM, x = {r,a,q} includes the spatial r and spin a variables 
and the index of nucleon type q = n,p. The integration over x includes the summation over a and 
Q- 

The functional depends on the total local density n{r). One can define other energy 
functionals which are dependent either on nq{r) or even on p{x,x). The particular choice of 
functional dependence is determined by the task under consideration. 

In the local DFT it is proved that the minimum of the functional E}{K[n] is just the true 
ground-state energy Egs and is attained for the true ground-state density ncs- To obtain Egs and 
ncS: one can use the KS method which yields single-particle equations similar to the HF equations. 
The rigorous derivation of these equations is based on the following statement proved by Lieb |14|: 

If n{r) > 0, / n{r)dr = N, J ^n{r))'^dr < oo, then there exists N -particle Slater- 
determinant wave function built up from an orthonormal set of N single-particle wave functions 

^o{xi,...,xn) = {Nl)-^/^det{ipi{xj)} , (7) 

such that ^0 ~^ n(r). In other words there is a many-to-one mapping of A^-particle Slater deter- 
minant wave functions onto the local particle density n(r). 

This theorem enables one to define the kinetic-energy functional To[n] for a system of 
non-interacting particles: 

To[n] = inf < ^olTl^-o >, (8) 
and to divide the HK functional EhkIit] into two parts: 

EHK[n] = To[n] + W[n\, (9) 

where the energy functional W[n] contains the potential energy as well as the correlation part of 
the kinetic energy. 

Because the density n is produced by some Slater-determinant wave function one has from 
Eqs. dD-®: 

N 

n,ir)=Y.J2\^,{r,a,q)\\ (10) 

i=l cr 

By the same reasoning the kinetic-energy functional of the non-interacting system (|8|) can be written 
as 

^oN=E^E / \^V^ir,a,q)\'dr. (11) 

i=l u,q 

Notice that one could define the kinetic-energy functional on the basis of a more general 
set of A^-particle wave functions: 

T[n] = inf < ^'|r|^' > . (12) 



However, this functional cannot be written in the form (|ll|) and thus it is useless to derive KS 
equations. 

Applying the variational principle to the functional EhkW] with functions (pi, if* as 
functional variables, one obtains in accordance with Eqs. @), (|9|)"(lTl|) the following KS equations: 

hHKfi = eifi , (13) 
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with 

hHK = -i^^ + U{r), (14) 
Zm 

where U{r) = 5W/6n is the local mean-field potential and Ei are the Lagrange multipliers to ensure 
the normalization condition of the single-particle wave functions (pi. 

Often the energy functional W[n] is divided into two parts: W[n] = Enin] + Exc['>T']i 
where EH[n] is the "direct" (Hartree) functional, while Exc[''t] is the exchange-correlation energy 
functional. Consequently, the mean field potential U is also divided into two parts. For the sake 
of simplicity we will not do this in the present paper. 



Eq. (14) does not contain a radial dependent effective mass nor a spin-orbit potential 
which are essential ingredients of the model nuclear single-particle Hamiltonian. To include them 
we propose the following method based on a special version of the non-local extension of the DFT. 
Let us define the energy functional: 

£o[po]= inf <^'o|^|1'o>, (15) 

where is any Slater-determinant wave function of the form (|7|), po is the single-particle DM 
produced by ^'o according to Eq. (^) (i.e. Slater-determinant DM) and H is an effective many- 
body Hamiltonian which generally does not coincide with the microscopic Hamiltonian H. In our 
approach the operator H plays the role of an arbitrary reference point, the choice of which will 
be discussed below. We have to note that at the present moment H is an arbitrary A^-particle 
operator such that the matrix element in ( p^ exists. 

The functional £^o[po] h^-s a form of the HE energy functional built up on the base of the 
effective Hamiltonian H. So in what follows we shall refer to it also as the HE energy functional. 
Let us define the residual correlation energy Erc as follows 

ERc[h]=E[h]-Eo[h], (16) 

where h = {up, n„} and 

E[h] = inf < ^H\^ >, (17) 

Eo[n] = inf < ^ol-H'l^o >= inf inf < ^'o|-H'|^'o >= inf £o[po] ■ (18) 

The quantity E[n\ is the exact energy functional built up with the true microscopic Hamiltonian 
([l|) on the set of any normalized wave functions ^ . The auxiliary functional E^lh] (as well as the 
kinetic-energy functional To[n] in the KS theory) is defined according to the Lieb theorem for any 
(not very "bad") local density n. The final energy functional of our version of the non-local DET 
is defined as: 

£[po\=£o[po] + ERc[n], (19) 
where po is related to h through Eqs.(§)-(|6|). The functionals <?o[po] and E^cin] are defined by 



Eqs. (15) and ([Tq). Eor the moment we shall not discuss if these functionals are known or not. The 
most important for us is that they are rigorously defined. 

The main property of the functional £[pq\ is expressed by the following equalities: 

inf £[po]= inf inf £ [p^] = inf E [h] = Eqs , (20) 

Po n po^n n 

where Eqs is the true ground state energy of the interacting system as in the case of the HK theory. 
To obtain the equations of motion we have to suppose that the choice of the Hamiltonian H in 



Eq. ( pISD ensures that the infimum of the functional £[pq\ in (20) is a minimum. In addition, we 
use the general formula for the Slater-determinant DM pQ under the variation which follows from 
Eqs. dD and (0). Namely 

N 



po{x,x') = Y,^i{x)^*{x'), (21) 
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where the sum is taken over the occupied states. Applying the variational principle to the functional 
S[po] defined by Eq. ( ]I9| ) and using the functions ipi, if* as functional variables according to Eq. (^) 
one will get the following set of single-particle equations: 

ho{x,x')ipi{x')dx' + URc{x)(pi{x) = eiipi{x), (22) 

where we have defined the non-local Hamiltonian /iq and the local potential Urc ^ follows: 

hoM = (23) 
dpo[x',x) 

Urc{x) = URc{r,q) = ^^^. (24) 

dnq{r) 

It is worthwhile noting that the occupation numbers of the Slater-determinant DM are fixed to be 
either 1 or 0. Thus, we avoid difficulties arising from the variation over the occupation numbers 



that appear in the theory developed in Ref. |13] 



3 REDUCTION TO THE QUASI-LOCAL THEORY 

The approach described above enables one to introduce a reduced energy functional which 
depends on the following set of local quantities: the local particle Uq, kinetic-energy Tq and spin Jq 
densities for neutrons and protons: 

nq{r) 



T„{r 



= ^ /"d2;'5(x-x')po(a;,x'), (25) 

= ^ I dx'8{x-x'){VrVr')po{x,x'), (26) 

Jq{r) = iY, j dx'5{r -r')5q^q,[{a)^,^^xV,]pQ{x,x'), (27) 

where 5{x — x') = 5{r — r')5a-^a'^q,q' ^ the quantities Tq and Jg are the uncorrelated neutron 
and proton kinetic-energy and spin densities respectively. Introducing the short notation pql = 
{rip, nn,Tp, Tn, J p, J n} , let US define the quasi-local energy functional as follows: 

[pql] = S^"- [pql] + Erc [h] , (28) 

where 

£0"^^= inf £o[po]. (29) 

PO^PQL 



Notice, that the many-to-one mapping po — > pql is established according to Eqs. (p^)-(|27|), and 
that the set n = {np,n„} enters Pql- n G Pql- 

From Eqs. (|o|), (|§ and (||) we have 



inff«^[pQi] =fe. (30) 

PQL 

Using Eq. (10) and the explicit expressions for the remaining local quantities: 

N 

rqir) = ^^|V(^,(r,a,g)|2, (31) 

i=l cr 

N 



1=1 , 
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one can apply the variational principle to the energy functional £^^[pql] with ipi, (p* as functional 
variables. The resulting single-particle equations are: 



hgipi = enpi , (33) 



where 

2m* (r) 



hg = -V——^V + Ug{r)-iWq{r)-[Vx(T], (34) 



and 



2m* (r) STq{r) 



(35) 



If the functional £'^^ were known, one would be able to calculate the exact ground state energy 
£^0 = Egs and exact local densities Hq = Uq^cs- The kinetic-energy density Tq and spin density Jq 
correspond to the system without correlations and do not coincide with the exact densities in our 
approach. 

The following remark is in order. As was pointed in the Introduction, the energy functional 
dependence on the exact kinetic-energy and spin densities could be introduced, for example, by the 
above procedure of the quasi-local reduction being applied to the energy functional of the non-local 
theory developed in Ref. [13|. However, the mapping of the Slater-determinant wave functions 



onto the exact kinetic-energy and spin densities is not defined, as opposed to the mapping onto the 
exact local density, i.e. Eqs. (|3l| ) and (^2|) are wrong for the exact densities. As a consequence, it 
is impossible to derive the single-particle equations such as (|33|) . The use of the Slater-determinant 
DM (po in our approach) enables us to avoid this problem. 



4 EXTENDED THOMAS-FERMI APPROXIMATION 
IN THE QUASI-LOCAL DFT 

Now we are turn to the discussion of the choice of the effective operator H. First, the energy 
functional £q has to be well defined. While the functional E[n\ is always defined, this is not the 
case for the functionals <fo[po] and -Eo[^]- Indeed, the matrix element of the microscopic Hamiltonian 
with the bare nucleon-nucleon (NN) interactions over the Slater-determinant wave functions can 
diverge due to the short-range singularity of the force as it is happens in the standard many-body 
theory. Thus, we will use a Hamiltonian H with an effective NN interaction whose matrix elements 
are well defined. Second, the local energy functional E[n\ obviously has the minimal property 
that is necessary to apply the variational principle; however, it can be not true for the non-local 
energy functional £ [po] (in spite of the fact that the equality (^) is always true) . Hence one has to 
choose H to ensure the minimal property of the energy functional ^[po]- Notice that at least in one 
particular case, when H = T, this condition is fulfilled because we come to the usual KS theory. 
Indeed, if one sets if = T then m* = m, Wq = 0, and Eq. (|33|) coincides with the KS equation. 
In this case the residual correlation energy functional Ejic[f>] corresponds to the sum of the direct 
(Hartree) term and the exchange and correlation energy in the KS method. 

Here we consider a more general case, setting H to be the A^-particle Hamiltonian with 
an effective NN interaction: 

fj = T + J2vf/ + J2vf/'^K (38) 
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We will use the effective NN force of the form 



^rj"" = vt,+v^, (39) 
where the central part of the effective force is given by 

Vt, = EK + bnP^j - hnP[j - mnPr,P[j]Vn{s), (40) 
n 

Wn,bn, hn,mn are the parameters of the force (n = 1,2,...), Pfj and P[j are the spin and isospin 
exchange operators and Vn{s) {s = ri — rj) are the radial formfactors of the central part of the 
effective force. For the sake of simplicity we will consider only one term in the sum (|^) i.e. we 
will drop index n in the following. The spin-orbit part of the force is chosen as follows 



vt° = iWoiai + <Tj) ■ [k' X 6{n - rj)k] , (41) 



where k = ^(Vj — Vj) denotes the operator acting on the right and k' = — ^(Vj — Vj) is the 
operator acting on the left. 

It was mentioned that the choice of the effective Hamiltonian H and the effective in- 
teractions is rather arbitrary. Practically the operator H is confined only by the above formal 
mathematical conditions. The situation is quite different in the usual HF method where there are 
no other ingredients apart from the effective forces which are taken to be density dependent in 
order to ensure the nuclear saturation. In our approach the particular choice of H is compensated 
by the addition of the residual correlation energy Efic[n], which contains all necessary density 
dependence of the total energy functional. In the applications of the method the functional Ejic[n\ 
is parametrized phenomenologically, the parameters are adjusted to describe nuclear ground-state 
properties. Following this ideology the effective interactions entering operator H are taken to be 
density independent. Thus, we would like to stress that the effective interactions in our DFT 
approach are not exactly the same as in HF theory. 

Let us define the density Tio of the quasi-local HF energy functional £q^[pql] as follows 

E^"- = j drHoir) . (42) 
According to (^)-(|4l|) the energy density T-Lq is described by six terms: 

2m ^ ^ ^^2:c?t + ^Dir + rLExch+ n- ■ (4cij 

The direct nuclear energy density H^l^^ comes from the central part of the NN force and is given 

by 

n'^fir) = \j dr'{{w + \)n{r)n{r') - {h + |)[n,(r)n,(rO + n„(r)n„(r')]}t;(|r -r'\). (44) 
The density of the Coulomb direct energy is 

Wg?."'(-) = y / (45) 

These direct energies give the contribution to the so-called Hartree functional. To calculate the 
exchange terms that come from the central part of the NN force we use the recently proposed 
ETF approximation for the DM up to 'h?' order [^. Notice that there are other possible options to 
obtain the quasi-local energy functional based on the Negele-Vautherin and Campi-Bouyssy DM 
expansions [Q, ^, ^ . In our approach for spin-saturated nuclei the nuclear exchange energy density 
is given by two terms 

IjNucl njNucl < njNucl I ACi\ 

i^Exch — i^Exchfl "T '^Exch,2 ■ V^OJ 
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The first term is calculated to order (which corresponds to the Slater approximation for the 
DM): 

n^TcU^) = J dsv{s) (^^Xei J2 (n,{r)UK^)f + X,2nn{r)h{Ks)np{r)Ukps)^ , (47) 

where kq{r) = (37r^ng(r))-'^/^ is the Fermi momentum, ji{x) = 3ji{x)/x, ji{x) is the spherical 
Bessel function and X^i = m + h/2 — h — w/2, Xe2 = m + h/2. The second term corresponds to 



the fi^ correction: 



WK,2(r) = E |^[(/. - - \kln, - ^An,) + _ JLa^,)]. (48) 
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In this equation fq = fq{r, kg), = {dfq{r, k) / dk)k=kq- The function fq{r, k) is the inverse of the 
position- and momentum-dependent effective mass and is defined in the ETF approximation by 

fir/:) II "' ^'^S.C-.^ 

where V^^Jj^^ is the Wigner transform of the exchange potential in the Thomas-Fermi approxima- 
tion: 

VExctpir, k)=J dse-'^^v{s)[X,inp{r)ji{kps) + Xe2n„(r)ji(A:„s)] (50) 

and analogously for V^^Jj^^ with the permutation of indices p and n (see |^ for details). It is 
worthwhile noting that within the semiclassical ETF approximation the kinetic-energy density is 
a functional of the local density. Then the energy functional obtained would only depend on the 
local particle density and spin density. However, it was found in |^] that the use of the quantal 
kinetic energy in Eq.(^), which yields the radial-dependent effective mass, improves agreement 
with results of the full HF calculation significantly. It motivates us to use the anzatz (^T]) for Tq in 
the present paper. 

The Coulomb exchange energy consists of the Slater term and the second-order correction 
that in the ETF approximation is written as [p^ , |8| : 

'HE.cdr)- ^^^J 432^(3^2)1/3 ^4/3 ' ^^^) 

Finally, the spin-orbit energy density is given by 



= _ Wo[n{r)VJ + n„(r)VJ„ + np(r)V Jp] , (52) 



where J = Jp + J n- 

It is worthwhile noting that in this section we replace the exact quasi-local functional 
by the approximate functional calculated within the ETF approximation. The difference between 
them gives a very small contribution (see the following sections), but it cannot be totally included 
within the residual correlation term because it depends on the pqi while the latter only depends 
on the h. The determination of the residual correlation energy ERc[n] will be specified in the next 
section. 

5 THE ENERGY FUNCTIONAL AND SINGLE-PARTICLE 
EQUATIONS FOR THE GOGNY FORCE 



The formulae (p^, (^6|)-(^0[) are valid for any radial formfactor v{s) of the central part of the 
effective forces. Here we give the explicit expressions for the above-defined quantities in the case 
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of a Gaussian formfactor v{s) = exp{—s'^/a'^) entering in the Gogny-type effective forces which are 
used in the numerical apphcations of our method. Assuming spherical symmetry of the particle 
densities we get from Eqs. and 



^Dt"'ir) = dr'r'iexp l-'-—^ - exp } 



h Tfi 

X + 2)"(^)™('^') - + y)K(^)"'p(^') + "'n(^)'^n(^')]} ' (53) 



37r5/2a3 



{Xei E ( ^a'k'^evfiak,) + ( ^ - 1 ) eM-a'k',) - + I 



q 



+ Xe2 E V[^a^{f^n + vkp){kl + kp-7]knkp)er{f^{kn+r]kp)\ 



.2 



+ (^a'^{kl + kl-r]knkp)-2^expi-^{kn+rikp)'^\]}. (54) 



The second-order correction to the exchange nuclear energy density (48) can be rewritten in the 
following way: 

where the explicit value of the functions Fg and Gq calculated with a Gaussian formfactor are given 
in Appendix 1. 



The SPP is defined following Eq.(36). According to ([2q), ( |4^ ) and (^3|) we split it into 
five pieces: 

TT TjNucl I TjNucl I TjCoul I TjCoul I ttRC (r.p.\ 

'-'q — '-^Dir,q + '-^ Exch,q + "^Dir.g + '-^ Exch,q + "^g ) l^OJ 

where the direct nuclear SPP is given by 

UDir%{r) = — (ir'r'{exp 1-^—^1 - exp 1-^—^1 } 

h 777 

X [^n, + -)n{r')-{h+-)nq{r')]. (57) 
The exchange nuclear potential consists of two parts following (^6|): 

jjNucl jjNucl I jjNucl 

^ Exch,q — ^ Exch,qfi ^ Exch,q,2 ■> \^°) 

where, for example, the Slater part of the exchange SPP acting on the protons is given by 

tNucI _ 2 r \- rV^„3,.3_r/ , A , „2,^2„,_/ „2,„2\ „2l2 



UexcKpA = -^=^{Xei [^a:'klexi{akp) + a'k^ exp{-a'k^) - aVp] 

+ 2Xe2 E ^[^«'A;3erf(^(fcp + r?A;„)) + ia2fc2exp(-^(fep + r?A;„)2)]}. (59) 



4 P ^ ' 2 P 4 

rj=±l 
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For the second-order contribution to the exchange SPP, we have 



UExch,p,2 = ^ ~ 5 P P> + 27 I 1|~ 





5 " 27 " Un 

^ (rP ^QFP^Ar^ ] ^ ^n(Vnp)(Vn^ 
+ 9FjAn.] - } 

-Fpk^p - ^AFp - ^AGp, (60) 

where the functions F^j' , G^', AFg and AGg used in Eq. (|60| ) are also given in Appendix 1. The 

formulae for the Slater and h'^ contributions to the nuclear exchange potential acting on neutrons 
are obtained by replacing n by p and p by n in Eqs.(^) and (|60|) . 

The Coulomb direct and exchange potentials, entering Eq. ([5^), are not equal to zero only 
for protons. In the explicit form we have 

t^2ctp(^) = (62) 

The inclusion of the /i^-correction to the Coulomb exchange energy (second term in Eq. (|5l|)) in 
the SPP leads to the unphysical behaviour of the potential at r — > oo. We thus only calculate its 
contribution to the binding energy as a perturbation. 

The radial dependent effective mass rn*{r) and the form factor Wq{r) of the spin-orbit 
potential are defined according to Eqs. (H), (|^, (pi), (E|), (|D and iM) by the relations 



2m* (r) " 2m^ ^^^^ 
Wgir) = ^Wo{Vn + Vng). (64) 

Now we shall turn to the determination of the residual correlation energy. We take it in 
the form of the phenomenological ansatz 

ERc[n] = jj drn^{r)[{2 + x^Wir) - (2x3 + 1) (n2(r) + nl{r))]. (65) 

The parameters t^, and a together with the parameters of the Hamiltonian H have to be 
chosen from the condition of the best description of the nuclear ground-state properties. In the 
calculations, which are presented in the next section, we use the well-known parameter set of the 
DIS Gogny force. Let us note that the formula ( |6^ ) is the standard ansatz which enters not only 
a density-dependent part of the Gogny forces but, for example, a density-dependent part of the 
Skyrme forces. Eq. ( |65| ) leads to the following contribution into the SPP in accordance with the 
definition (|56|): 



jjRC ^ ^ a-l 



(2 + a)(2 + X3)n^ - (2x3 + 1) [a{nl + n^) + 2ngn)] . (66) 
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6 NUMERICAL RESULTS 



In this section we want to check our local approximation to the HF method using finite-range 
Gogny forces. First of all, we compare the exact HF ground-state binding energies as well as the 
rms radii of the neutron and proton densities of some magic nuclei computed with the Gogny DIS 
force [18 1 with our DFT results. In this comparison we use two different quasi-local functionals : 
DFT-K\ where the exchange energy coming from the finite-range part of the interaction is taken 
at pure Thomas-Fermi level (Slater approach) and DFT-?l^ where the ETF-h? contributions have 
been added to the Slater part. Notice that in this DFT-^^ approach the semiclassical kinetic-energy 
density entering Eq. (48) has been replaced by the corresponding quantal density for the reasons 
pointed out above. In both DFT calculations we solve the local Schrodinger equation (^) for 
neutrons and protons with the potentials and effective masses reported in the previous section. 
Table |^ collects all these binding energies and radii which have been computed taking into account 
the two-body centre-of-mass correction. In our calculation we take into account this correction as 
explained in Appendix 2. We would like to say in passing that the numerical value of this two-body 
centre-of-mass correction along the whole periodic table is very well reproduced by using the pocket 
formula based on the harmonic oscillator and derived in |1S]. 

From Table [l| we can see that the DFT-?i^ binding energies reproduce the HF values fairly 
well. The differences between HF and DFT-^^ are smaller than 1% for all the considered nuclei 
from ^''Ca to ^''^Pb and in the case of ^^O the relative difference is only 1.8%. The DFT-h^ binding 
energies show larger discrepancies with the full HF results. The relative differences range from 7% 
in ^^O to 1% in the heaviest nucleus considered ^"^Pb. As regards the rms radii of the neutron 
and proton densities, the full HF values are again better reproduced by the DFT-?i^ approximation 
than by the DFT-?i'^ approach. These results show the importance of the corrections in the 
local approximation to the HF exchange energy. It should be pointed out that the eigenvalues Si 
in Eq. ( |3^ ) have no rigorous physical sense in the DFT except for the energy of the last filled level, 
which corresponds to the neutron or proton separation energy (chemical potential). Table || also 
displays the neutron and proton chemical potentials obtained using the DFT-;i°, DFT-^^ and HF 
approximations. The DFT-?i^ chemical potentials differ from the HF chemical potentials by less 
than 1 MeV while the shift of the T>FT-h^ separation energies with respect to the the full HF values 
is larger and can be around 3 MeV for light nuclei. 

It is important to note that the agreement of the proposed DFT approximations with the 
full HF results is determined by the treatment of the non-local effects. The contribution of the 
pointed effects can be quantified in terms of the effective mass in nuclear matter. The results of 
Table H show that for effective forces with an effective mass over the bare mass of around 0.7, as in 
the case of the Gogny forces [20|, the non-local effects can be very well accounted for the DFT-7i^ 
functional proposed in this paper. However, for forces where the non-local effects are larger, the 
result of our DFT-?i'^ approximation is worse when compared with the HF results, but is still better 
than other suitable choices of the exchange-energy localization such as the Negele-Vautherin or 
Campi-Bouyssy approaches (see Ref.[^ for more details about this point). 

Figures |l] and ^ displays the neutron and proton densities for ^''Ca and 208p]-, obtained 
using the DIS force in the full HF (solid line) as well as in the DFT-^^ (dashed line) and BFT-h° 
(short-dashed line) calculations. The DFT proton and neutron densities nicely reproduce the 
surface and the tail of the full HF densities. In the region near the centre of the nuclei, the DFT 
density follows the full HF density profile very well although a small shift between the full HF 
and DFT proton and neutron densities appears in this central region. These differences can be 
attributed to the fact that our DFT description of the s-orbitals, whose wave functions mainly 
provide the proton and neutron densities at the centre of the nuclei, show some small differences 
with the corresponding HF s-orbitals. Comparing the DFT-^'^ and DFT-?i^ densities, it can be 
seen that including the contributions in our local approximation, one obtains a better agreement 
with the full HF densities. 

Figures 3 and 4 display the radial dependence of the neutron and proton effective masses 
over the bare nucleon mass calculated with the T)FT-h^ approach (see Eq. (|63D) for ^''Ca and 
208p|-j j^ucigi (solid lines). Because there is no explicit radial-dependent effective mass in the full 
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Figure 1: Neutron and proton densities of the nucleus ''Ca calculated with the DIS Gogny force 
using the DFT-^° (short-dashed line) and DFT-^^ (dashed line) approaches compared with the full 
HF densities (solid line). 
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Figure 3: Neutron and proton radial dependence of the effective mass of the nucleus Ca calculated 
with the DIS Gogny force using the DFT-^^ approach (solid line) compared with the corresponding 
HF effective masses obtained with the Skyrme SIII force (dashed line). 



15 




2 4 6 8 10 12 

r (fm) 




2 4 6 8 10 12 

r (fm) 



Figure 4: Same as Figure |3| for the nucleus ^"^Pb. 
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Table 1: Total binding energies B (in MeV), neutron r„ and proton rp rms radii (in fm) and 
separation energies of neutrons 5„ and protons Sp (in MeV) of some magic nuclei computed with 
the DIS Gogny force using the DFT-h^ and DFT-h^ approaches compared with the full HF results. 
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HF calculation of finite nuclei using the Gogny forces, we compare the DFT-?!^ results with the 
neutron and proton effective mass over the bare nucleon mass obtained using the Skyrme III force 
|]2l|| (dashed lines). We find that the DFT-Tl^ results exhibit similar trends to those of the Skyrme 
effective masses. The differences between the two calculations are basically due to the different 
values of the nucleon effective mass in nuclear matter which are m* /m = 0.70 for the Gogny DIS 
force and m* /m = 0.76 for the Skyrme III interaction. 



7 SUMMARY 

In the present paper we propose a non-local extension of the DFT and its quasi-local reduction. 
To this aim we define an energy functional which depends on the Slater-determinant DM with the 
occupation numbers being either 1 or 0. This enables us to avoid the difficulties of the non-local 



DFT reported in [13|. Defining the uncorrelated kinetic-energy densities and spin densities we 
construct the quasi-local energy functional and rigorously derive the single particle equations with 
the radial-dependent effective mass and the spin-orbit potential. 

In order to define the energy functional of the Slater-determinant DM one has to introduce 
an effective Hamiltonian which ensures its existence. This feature of nucleon systems arises from 
the specific properties of the bare NN force in contrast to the Coulomb force in electron systems. In 
our approach the total energy functional consists of HF part and residual correlation energy. The 
HF energy functional can be calculated directly, while the residual correlation energy is considered 
phenomenologically. Using the recently proposed semiclassical ETF approximation for the DM Q 
we obtain a quasi-local energy density functional only depending on the local particle, kinetic-energy 
and spin densities. The resulting single-particle equations of motion contain the local mean- field 
potential, the uncorrelated effective mass and spin-orbit potential. Using the finite range density- 
dependent Gogny force they are calculated analytically. The use of a different effective force such 
as M3Y will be reported in a forthcoming publication. 
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We apply our method to the calculations of the nuclear ground-state properties and com- 
pare it with the exact HF solutions using the Gogny DIS force. We obtain a very good agreement 
in the description of the binding energies and root mean square radii. The single-particle energies 
of the highest occupied neutron and proton levels in the full HF calculation are well reproduced 
by our local approximation. The particle densities are also in good agreement with the exact HF 
densities. Wc analyze the radial dependent effective mass within our approach. Comparing it with 
the result obtained with the Skyrme HI interaction, it also exhibits a very reasonable behavior. 

In conclusion, our approach has the following advantages: it handles local differential 
equations in contrast to the intcgro-differential equations in the HF approach, at the same time the 
quality of the obtained results is sufficiently high; our method enables one to construct a quasi-local 
energy density functional on the base of the effective forces with arbitrary radial formfactors; the 
method can be straightforwardly generalized to the non-spherical case. 
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Appendix 1 



In this appendix we present the exphcit expressions for the functions Fq, Gq, , , AFq and 
AGq calculated with a Gaussian formfactor. These functions are used to obtain the second-order 
contributions to the exchange nuclear energy density ( ]55| ) and the corresponding SPP (pOl). In the 
following, excepting Eqs. ( [70| ) and (^), we assume that q' / q. 

F, = -^{Xeiz^xp(-|)Qi(|)+Xe24exp(-ili)Qi(^)}, (67) 

2 2 2 2 

Gq = ^zl{X,,zleM-'-^)[QlC-^)-4Q2C-^)] 

+ Xe24exp(-i±i)[Qi(^)-z2Q2(^)]}, (68) 
where Zq = akq and the functions Qm{x) are defined by 

sinh(x) 1 dQm{x) . . 

Qo{x) = — - — , ^ 2x — dx — ■ 

The functions F^' , G^', F^''^" and G^'^" (the two last functions needed to obtain AFq and AGq, 
see below) are defined as: 

pg' = J_^ pi'i" = ^^^1 (70) 

kq>dkq>' 5 dkq,dkq,r ^'^> 

, ^ _]_dGq ^q'q" ^ _d^Gq_ 

" V ^ dkq,dkq. ' 

and their explicit form is: 
40F 



F^ = --^{2XeieqZq[i3 - zl)Qi{xq) + 4^2(2:9)] + Xe2eo4 [z,2,Q2(a;o) - Qi{xq)]}, (72) 



= -^Xe2eov[(6 - 4)Qi(^o) + ^|4Q2(xo)], (73) 



a 



4 



= --^{4.XeieqZq[{2 - z„^)(3 - 2z„^)gi(xg) - 4(1 + 2z„^)Q2(x. 



Xe2eo4[(2 - - 4)Qi(xo) + 2^(4 + Zg^)Q2(xo)]}, (74) 
= -^Xe2eov[(24 - 144 + zlzl + 4)Qi(xo) + 2^4(2 - 4)^2(^0)], (75) 



pq'q' " 



^4 

= __Xe2eoZq4[(6 - 2z2,)Qi(xo) + 4(^2 + 4)Q2(2;o)], (76) 



Gl = ^{2Xeieg4[(5 - 2zl)Q^{xq) + 4(3 + 24)Q2(x,)] 

+Xe2eo4[(4 - 4 - zl,)Qi{xo) + 24,(3 + z2)Q2(2;o)]}, (77) 
q4 

= ^Xe2eo4v[(6 - 24)^1(2:0) + 4(4 + 4)Q2(xo)], (78) 
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+Xe2eo4[(8 - lOz^ + 44 + 3^24 + 4)Qi(xo) 

-4(48 + 10^2 + + 34)Q2(xo)]}, (79) 



+4(4^2 + 24 - 3^24 _ zi)Q2{xo)], (80) 



= ^^Xe2eoZy4 [(24 - 6zg2 - 84 + 3^24 + 4)Qi(xo) 
+4 (4^2 - 64 - 3^24 - 4)Q2 (xo)] , (81) 
where Zg = akq, Xq = Zg/2, xq = ZpZn/2, = exp(— Xq), cq = exp(— (xp + Xn)/2). Finally AFq and 



AGq are given by: 

7r2 1 CVn 12 

+-^[(F/^' - 2F/)^^^ + 3F/An,,] + ^i^f (Vn,)(Vn,0}, (82) 
tP- 1 fVr? 12 

= Y^fc^t^^'' " ^^^^^^ + ^^'^""^l 

+i-[(G^^' - 2G^')^^^^ + 3G^ An,,] + (Vn,)(Vn,0}. (83) 
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Appendix 2 

In this Appendix we describe briefly the method to calculate the centre-of-mass correction to the 
ground-state energy. As is well known, the general idea consists of subtracting the quantity 

p2 

(84) 

from Egs- Here P is the total momentum operator, M is the total mass of a nucleus. Usually the 
quantity E'-^'^ is represented as a sum of two terms: 

j^CM ^ j^CM ^ ^CM ^ (85) 

where is the one-body, £'2'^^ the two-body centre-of-mass kinetic energy. The quantity 

E^^ is defined by formulae: 

E?'' = E i^iT , i^iT = ^Sp (pV,) , (86) 

where in accordance with the definition (P) the following notation is introduced 

Pq = Pq{T^, r') = E q, r', a, q) . (87) 

(7 

Hereinafter the symbol Sp denotes the trace over the space variables. The subtraction of E^^'^ leads 
to the simple renormalization of the nucleon mass in the single-particle Hamiltonian hq-. niq fhi^g, 

fhi^g/mg = M/{M - nig) . (88) 

The reasonable method for the evaluation of the quantity E2^^ is the Hartree-Fock approximation 
for the ground-state wave function "^gs in Eq. ([8^). In addition we adopt the following approxi- 
mation for the single-particle DM: 



p(r,o-,g, r',a',q) = ^S^y Pqir,r') . 



With these assumptions we have 



E2 



CM _ 

q 



E E^,r ' = -^SP {K^fPq) , (90) 



where the single-particle operator K2^^^ is defined as 

K^f=^PPqP. (91) 

In contrast to the one-body contribution, the subtraction of £"2"*^ leads to the additional non- 
locality in the hg because in the self-consistent approach we have to add the non-local operator 
^2,q^ to the single-particle kinetic-energy operator. So the total correction to the hg is: 

■ +K^f. (92) 



2mg 2mi^g ^''^ 

In the local or quasi-local DFT and in similar approaches the non-locality of K2g^ in Eq. (| 



leads to unpleasant difficulties. So we use the simplified method, proposed in Ref. [16|, to take into 
account the contribution of the operator in the single-particle equations. 
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Let us write the density matrix pq in the form: 



/r/k 



(93) 



where R = (r + r')/2, s = r' — r, ng(R, k) is the momentum distribution function. The approxi- 
mation consists of replacing function ng(R, k) in Eq. (^) by some effective constant value . 



Substituting Eq. (^) with nq(R, k) 



n 



CM 



into Eq. (||) we get 



f?CM 



n 



cmP 

M 



(94) 



The value of n^^ is defined by the substitution of K^^^ into Eq. (^ instead of K^^^ . Taking 
into account Eqs. 



and (94) we obtain: 

:CM 



n„ 



CM I TpCM 



jpUM I TP' 



(95) 



The quantities E^^^ and E^^^ are defined everywhere by Eqs. (|86|) , (90) which can be rewritten 
in the following forms making use of Eq. (^) for the Slater-determinant density matrix 



CM 



-2 N 



T7"~ 



E. 



CM 
2,9 



i=l o" 
f^2 N N 

~Tm E E I E / dv'P*^^^ ^' 9)V(/Ji'(r, cr, q)\^ . 

i=l i' = l o- 



(96) 
(97) 



It is easy to prove, using Eqs. (|9q) , (|97D and the completeness of the set of functions (pi, that 
\^2,q^\ < l-^i^g^l aiid consequently: 

(98) 



One can also prove that in the limit cases (n?**^ 



and n^*'^ 



1) the action of the approximate 



operator K2^^ upon the wave functions of the occupied orbitals ipi coincides with the action of the 

initial operator Ki^^^ as defined by Eq. (91). So Eq. (94) can be considered as the interpolation 
formula. 

The total centre-of-mass correction to the hq in this method is reduced to the renormal- 
ization of the nucleon mass as in the one-body case: 



2mn 



2ma 



rrin 



M 



M + {2n^^ - l)mq ' 



(99) 
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